Non-classicality of photon added coherent and thermal radiations 
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Production and analysis of non-Gaussian radiation fields has evinced a lot of attention recently. 
Simplest way of generating such non-Gaussians is through adding (subtracting) photons to Gaussian 
fields. Interestingly, when photons are added to classical Gaussian fields, the resulting states exhibit 
non-classicahty. Two important classical Gaussian radiation fields are coherent and thermal states. 
Here, we study the non-classical features of such states when photons are added to them. Non- 
classicality of these states shows up in the negativity of the Wigner function. We also work out 
the entanglement potential, a recently proposed measure of non-classicality for these states. Our 
analysis reveals that photon added coherent states are non-classical for all seed beam intensities; 
their non-classicality increases with the addition of more number of photons. Thermal state exhibits 
non-classicality at all temperatures, when a photon is added; lower the temperature, higher is their 
non-classicality. 



I. INTRODUCTION 



Several branches of quantum optics from non-linear 
optics to laser-physics and cavity QED are very actively 
engaged in a variety of processes producing non-classical 
light. Such radiation fields attract attention, not only 
because they provide a platform for testing fundamental 
concepts of quantum theory, but also for applications of 
importance like precision measurements in interferome- 
try Q- Moreover, rapidly developing area of quantum 
computation and information theory has kindled further 
interest in generating and manipulating non-classical 
radiation fields - called quantum continuous variable 
states. These states are promising candidates for marry 
applications of quantum information technology 0, 0. 
In such a context, Gaussian light fields gain prominence, 
both in view of their conceptual and experimental im- 
portance. However, a need to leap beyond Gaussian do- 
main has been emphasized Q| and the degaussification 
process has been catching a lot of interest. Degaussi- 
fication can be realized in a simple manner by adding 
(subtracting) photons to (from) a Gaussian field and 
the resulting states are known to exhibit non-classical 
properties such as negativity of the Wigner function [j| , 
antibunching 6], sub-poissonian photon statistics or 
squeezing in one of the quadratures of the field Q etc. 

Almost a decade ago Agarwal and Tara introduced, 
theoretically, a new class of non-Gaussian states, which 
is obtained by repeated application of the photon cre- 
ation operator on the coherent state. The resulting class 
of states were identified to lie between the Fock state 
and the coherent state and were indeed non-classical. 
Recently, single photon excitation of a classical coher- 
ent field has been generated experimentally [To| and 
ultrafast, time-domain, quantum homodyne tomogra- 
phy technique has explicitly demonstrated a quantum to 
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classical transition. In another development, a traveling 
non-Gaussian field was experimentally produced by sub- 
tracting a photon from a squeezed vacuum £j. While 
the pulsed homodyne detection scheme confirmed non- 
Gaussian statistics for the photon subtracted squeezed 
vacuum, the Wigner function reconstructed from the 
experimental data failed to exhibit negativity. Kim 
et al. |ll| analyzed the non-classicality of photon sub- 
tracted Gaussian fields and identified that the photo de- 
tection efficiency as well as the modal purity parameter 
were not high enough to record a negative Wigner func- 
tion in the experiment of Wenger et al. Moreover, 
Kim et al. show that unless the input Gaussian radiation 
is non-classical, one cannot generate a non-classical field 
through photon subtraction 0. They also point out 
the contrasting situation of photon addition to Gaus- 
sian fields, where even a highly classical state like ther- 
mal state turns out to be non-classical [l2L Il3| . It may 
be noted that photon addition to a thermal state results 
in the removal of the vacuum part 01 an d hence leads 
to a truncated thermal state. All such truncated states, 
with their vacuum contribution removed, are shown to 
be non-classical by Lee ^3 • 

In view of the current experimental progress 0, IToj 
in the production of non-Gaussian radiation fields, and 
also in the tomographic reconstruction 0, 0, Q of 
Wigner functions of quantum states, it is timely to ana- 
lyze the non-classicality of photon added classical radi- 
ation fields, through negativity of their Wigner func- 
tions. In this paper, we investigate photon added coher- 
ent and thermal states. 

Entanglement, another striking quantum feature, has 
occupied a central position in the development of quan- 
tum information processing 0, 0- There has been a 
considerable progress in understanding the connection 
between non-classicality and entanglement. It has been 
identified 0, ^| that non-classicality is an unmistak- 
able source of entanglement. A beam splitter is capa- 
ble of converting non-classicality of a single mode radi- 
ation into bipartite entanglement. This property viz., 
non-classicality as an entanglement resource, has been 
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employed recently |l5j , to identify Entanglement poten- 
tial (EP) - a computable measure of non-classicality - 
of single mode radiation fields. EP allows us to ana- 
lyze the degree of non-classicality of a given single mode 
radiation field. We compute EP of photon added co- 
herent states (PACS) and show that the EP reduces 
with the increase of seed beam intensity, which is in 
confirmation with the analysis of the Wigner function 
of the state. A comparison of EP's of single, two and 
three photon added coherent states reveals that non- 
classicality of PACS increases with the addition of more 
number of photons. We verify that the Wigner func- 
tion of a photon added thermal state is negative at the 
phase space origin for all temperatures, while the EP of 
the state, evaluated in the low temperature limit, reveals 
that non-classicality of the state reduces with increasing 
temperature. 



II. 



MEASURES OF NON-CLASSICALITY 



Generally a non-classical state is recognized as one, 
which cannot be written as a statistical mixture of 
coherent states. It has been well accepted that the 
non-existence of a well defined Glauber-Sudarshan P- 
function ^3 implies non-classicality of a given state. 
However this identification poses operational difficul- 
ties, as it requires complete information of the state 
to be examined, so that it's P-function can be recon- 
structed. Several operational criteria, which are equiv- 
alent to the one based on the P-function and which 
can be used to distinguish between classical and non- 
classical states in experimental measurements have been 
proposed from the early days of quantum optics. Such 
signatures of non-classicality, verifiable in a simple ex- 
periment are, antibunching la and sub-poissonian pho- 
ton statistics 0, squeezing photon number oscilla- 
tions 0| , negative value of Wigner function H , etc. 

Here, we focus our attention on the Wigner function of 
a given quantum state. The non-classicality character of 
a state is strongly registered by negativity of the Wigner 
function. Especially, Fock states show a negative dip 
around the phase space origin, as has been clearly re- 
flected in the experimentally reconstructed Wigner func- 
tion [3- Moreover, there has been an ongoing effort 
towards more efficient quantum homodyne tomographic 
techniques P. ITolll^| . and in such a context, analysis of 
Wigner function proves to be useful. 

The Wigner function of a system, characterized by 
the density operator p is defined [{| through 

W (^p) = \j^ + y\p\i - y) e ~ 2tpy d y ■ 0) 

Basically Wigner function is a quasi-probability distri- 
bution representing quantum states in phase space. It 
is not a true probability distribution as it can take neg- 
ative values also. If, for a state, Wigner function takes 



negative value, the quantum state has no classical ana- 
log. However, the converse does not hold good: When 
the Wigner function is positive everywhere, one can not 
conclude that the state is classical. For example, for 
a squeezed state, Wigner function is a Gaussian and is 
positive throughout. But, squeezed radiation is one 
of the most important non-classical field. Thus, one has 
to resort to other measures of non-classicality. 

There have been several approaches to quantify non- 
classicality of a state through universal measures like, 
Hillery's non-classical distance [l9( and Lee's non- 
classical depth [2(j. Distance between a given non- 
classical state and the set of all classical states is non- 
zero and hence, serves as a measure of non-classicality, 
called non-classical distance |19|. However, identifying 
an optimal reference classical state is one of the main 
problems associated with these distance based mea- 
sures [2ll| of non-classicality. 

Lee's |20j non-classical depth, < r m < 1, is es- 
sentially the quantity of smoothing required to trans- 
form a non-positive P-function into a well behaved pos- 
itive distribution. (Non-classical depth r m is also in- 
terpreted as the minimum average number of thermal 
photons that are necessary to destroy the non-classical 
effects of a given state). A classical state has r m = 0. 
For a pure Gaussian state (Squeezed vacuum), the non- 
classical depth varies between and \ [3> for all non- 
Gaussian pure states r m = 1 jU, thus placing all such 
states to be identical, as far as their non-classicality 
is concerned. Therefore, this criterion forces one to 
conclude that the non-classicality of PACS is indepen- 
dent of the seed beam intensity. Moreover, non-classical 
depth of photon added thermal states (for that mat- 
ter, all truncated states) is shown to be a maximum i.e., 
T m == 1 E3- According to this measure, photon added 
thermal states are equally (maximally) non-classical at 
all temperatures. 

In this paper, we consider Entanglement potential 
(EP) [3> as a universal measure of non-classicality 
for our discussion of PACS and photon added thermal 
states. EP gives the amount of two-mode entanglement 
that can be generated from a non-classical input state, 
in a linear optics set up. It is important to note here 
that a classical single mode radiation does not get en- 
tangled in such an arrangement [l6l | . EP is nothing but 
the logarithmic negativity |23| of a bipartite quantum 
state Qo , which results from mixing a given single mode 
state tr, with vacuum state, in a 50 : 50 beam splitter. 
More specifically, EP is defined as 



EP = log 2 ||^ T Hi, 



(2) 



where g^ T denotes the partial transpose of a two-mode 
density operator p a = Ubs{°~ ® \fy{Q\W B s- In equa 



tion (J2J, the symbol || • ||i denotes the trace norm[2J 
and Ubs corresponds to a 50 : 50 beam splitter i.e., 
Ubs = exp(|(a t 6 — atf)), whose action (Heisenberg 
view point) on the creation operators a) ,w of the two 



input ports is explicitly given by 

Ubs JUbs = -^(o t + 6 t ) 
U BS & 1 U BS = (6* - 



(3) 



TTie state cr is said to &e non- classical, iff its entangle- 
ment potential is nonzero. EP has been evaluated [l5j 
for a variety of non-classical states like squeezed states, 
even and odd coherent states, Fock states etc. 

The remainder of this paper is devoted to a study 
on photon addition to (i) coherent state, an example of 
a pure Gaussian state and (ii) thermal state, a mixed 
Gaussian state, both of which are well-known classi- 
cal states. Action of photon creation operator on these 
states results in non-classical, non-Gaussian states. We 
study the non-classicality of these states through nega- 
tivity of the Wigner function and the entanglement po- 
tential. 



Photon added coherent state 





Coherent states are the analogs of classical radiation 
fields. These states are described by a Poissonian photon 
number distribution and have a well defined amplitude 
and phase. It is interesting to see how these states turn 
non-classical, when a single quantum of radiation excites 
them. 

Recently ^(J; single photon added coherent states 
[SPACS] has been generated experimentally and to- 
mographically reconstructed Wigner function for such 
states has been analyzed. SPACS are obtained by ap- 
plication of creation operator a* on a coherent state \a). 
Normalized SPACS is given by 



| SPACS) = 



a^\a) = a^Dq{a)\0a) 

(i + M 2 ) V(! + H 2 ) 

= D a (a)[\l a )+a*\O a }} 

v/(i + M 2 ) 



, (4) 



where D a (a) = exp (a^a — aa*) is the displacement op- 
erator and the coherent state |a) = D a (a)\O a ); |0 a ) de- 
notes the vacuum state. Here we have used the property, 
Dl(a) aW a (a) = a^+a*. Wigner function of a SPACS 
is given by 



W(q,p) 



-£i(|2e-a| 2 ) 
^i(-M 2 ) 



exp(-2|c- a | 2 ), (5) 



where c = -^=(q + ip) and L\{z) = 1 — z is Laguerre 
polynomial of first order. Note that W(q,p) of the 
SPACS is negative when \2c - a\ 2 < 1. 

In Fig. ^ we have plotted the Wigner function of 
SPACS, for various beam intensities \a\ 2 . It is clear from 
these plots that negativity of Wigner function reduces 
with the increase of the intensity \a 2 \. 




FIG. 1: Wigner function of SPACS for different beam inten- 
sities \a\ 2 . For simplicity, we have chosen a to be real. Here, 
(a) q = 0.1 (b) for a = 0.9 and (c) for a — 3. 



We may verify whether addition of more number of 
photons leads to higher non-classicality of the beam. To 
see this we consider two photon added coherent states. 
Wigner function of such states is explicitly given by, 



W(q,p) 



L 2 {\2c- 



7rL 2 (-\a\ 2 ) 



exp(-2|c-a| 2 ), (6) 



2z 



is Laguerre polynomial of 



where L 2 (z) = 1 
second order. 

Fig. gives the plots of the Wigner functions of two 
photon added coherent state for varying seed beam in- 
tensities. We note the same behavior here too, viz., the 
non-classicalities - depicted through the Wigner func- 
tions - decrease with increasing intensity \a\ 2 . But 
through these plots we can not conclude if addition of 
more number of photons leads to higher non-classicality 
or not. 

To investigate this, we now evaluate the EP of these 
photon added states. Let us first consider a SPACS. 
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sity operator is given by 

A = WW 



with 



go 



1 



2(1 + |«| 2 ) 



V 



(9) 



/ 2 |a| 2 ^a* V2a* • • \ 

V2a 1 1 

V2a 1 1 





(10) 



in the Fock state basis \n a mj); n a , vtih — 0, 1,2, • • 

Since g' Q is locally equivalent to go (as D a -D& (^75 

corresponds to a local displacements on the two mode 
states), EP of g' Q and that of go are same. So we proceed 
with the evaluation of EP of the state go itself. 
The partial transpose of go is given by 



»(q.p) 0.2 




FIG. 2: Wigner function of two PACS for different beam 
intensities \a\ 2 . For simplicity, we have chosen a to be real. 
Here, (a) a = 0.1 (b) for a = 0.9 and (c) for a = 3. 



When SPACS is mixed with vacuum state |0) and is 
sent through a 50 : 50 beam splitter, the resulting two 
mode state is given by 



= [/ BS (|SPACS> <g> |0» 



V2j 



|l a 0|,) + |0al6)+V2a*|0a0i 

V/2U + H 2 ) 



since a 50 : 50 beam splitter Ubs acts on D a (a) as (see 
equation 

C/S5 A,(C0 4s ^ A, (^) D " (^) > ( 8 ) 

with A>(^=) = exp(^(at a -aa*)) and £> 6 (^) = 
exp (A=(b'a — ba*)). The corresponding two mode den- 



Q P o T 



2(1 



at 



/2|a| 2 


V2a 


V2a* 
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• 




v^a* 


1 








• 




V^a 





1 





• 
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Q 
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■ 




v • 










■J 



(11) 

where, except for the first 4x4 diagonal block all the 
other elements are zero. The non-zero eigenvalues of 
gg T are easily identified to be, 



A a = 



^3,4 — 



1 



2(1 



A 2 = - 



1 



2(i + M 2 r 



(i + | Q | 2 )± v /(i + | Q p)^-i) 

2(l + |a| 2 ) 



(12) 



Note that A2 is the only negative eigenvalue of g PT , and 
hence, the EP of a SPACS is given by, 



(7) 



EP = log 2 ||^ i || 1 =log 2 (l-|-2|A 2 |) 



= l0g2 



1 



a 



(13) 



Following similar lines, we can evaluate entanglement 
potentials of two and three photon added coherent states 
also. But the expressions are lengthy and do not ex- 
hibit a simple structure. We have computed them nu- 
merically and plots of entanglement potential for single, 
two and three photon added coherent states are given 
in Fig. |3| From the figure it is evident that EP is non 
zero for low intensity of the seed beam and it reduces 
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FIG. 3: Entanglement potential for single, two and three 
photon added coherent states. (Here, a is chosen to be real 
for simplicity) 



gradually with the increase of intensity showing that 
the state is non-classical for all seed beam intensities. 
This observation is in confirmation with the conclusions 
reached through Wigner function analysis. While the 
EPs of single, two and three photon added coherent 
states converge for higher beam intensities, they are all 
different for low beam intensity, with larger value for 
higher photon added states. This observation reveals 
that non-classicality of photon added coherent state in- 
creases with the addition of larger number of photons. 



B. Photon added thermal state 

Density matrix of single mode thermal state of system 
in thermal equilibrium, characterized by the Hamilto- 
nian H — a* a hw is given by 



10 -10 



exp 



Pth 



(-*) 



Tr 



exp 



(-*) 



(14) 



In the Fock state basis, pth can be expressed in the form 



Pth = A Y,x n \n)(r 



(15) 



where A = 1— x, x = e~~ ; < x < 1. Note that x — > 
limit corresponds to T — > and x — > 1 implies T — > 
oo. Photon added thermal state is obtained through the 
application of creation operator on the thermal state i.e., 



a f Pth a 



= (l-x)J2x n a^\n)(n\a 
n=a 

oo 

= (l-x)Y / X n (n+l)\n+l)(n + l\ 



n=0 



(16) 




FIG. 4: Wigner functions of single photon added thermal 
state for various temperatures: (a) x = 0.1 (b) x = 0.5 
(c) x = 0.9. 



Simplifying equation Ijlril) we get, 



\m}(m\ 



(l-x) 



d 



i 



■ pth 



dx \ 1 — x 

Normalized photon added thermal state is given by 

d ( 1 



(17) 



pa 



Pth 



dx \ 1 — x 



Pth 



(18) 



Wigner function of a thermal state has been identified 
to be 



W t h(q,p) 



-B exp [-B(q 2 +p 2 )}; B = |— -. 

7T 1 + X 



(19) 



Making use of the above equation, it is easy to evaluate 
the Wigner function of a photon added thermal state: 



1 



dx \ 1 — x 



Wth(q,p) 



(20) 



After simplification, we get the Wigner function of the 
photon added thermal state as 



W^(q,p)^-B 2 

TT 



2(q 2 +p 2 ) 



L exphB(q 2 +p 2 )]. 

(21) 

It is clear that the Wigner function is negative 

at the origin of phase space, at all temperatures. We 
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have plotted the Wigner function for different values of 
the parameter x - which in turn corresponds to various 
temperatures - in Fig. It is clear from the plots that 
the photon added thermal states at various temperature 
are all non-classical. 

Entanglement potential of photon added thermal 
states can be evaluated in the low temperature limit, 
since the higher Fock states have lesser occupancy in 
this limit leading to the truncation of the Hilbert space. 
Retaining terms up to first and second order in the pa- 
rameter x, the density matrices of photon added thermal 
state are given below: 

I order in x : ff£ = (1 - 2x) + 2x |2)(2| 

II order in x : ^ = _L_ [(1 - 2a:)|l)<l| 

+2x(l - 2x)|2)(2| +3x 2 13)(3|] 

(22) 

Following similar steps to evaluate EP of SPACS, we 



can numerically compute the entanglement potentials 
of photon added thermal states too, in the low tempera- 
ture limit. With the help of this analysis we realize that 
EP reduces with increasing temperature. We therefore 
conclude that non-classicality of photon added thermal 
states reduces gradually with the increase of tempera- 
ture. 

In conclusion, we have analyzed the non-classicality 
of photon added coherent states and thermal states, us- 
ing (i) negativity of the Wigner function and (ii) entan- 
glement potential. We have shown that photon added 
coherent states are non-classical for all seed beam in- 
tensities; the degree of non-classicality reduces with the 
increase of intensity. Photon added thermal states are 
shown to be non-classical at all temperatures and their 
non-classicality reduces with the increase of tempera- 
ture. 

We thank J. K. Asboth for useful comments and the 
Referees for their insightful suggestions in the light of 
which the paper has been revised. 



[1] Quantum Interferometry, edited by F. De Martini et al. 

(VCH, Weinheim, 1996). 
[2] W. Tittel, G. Ribordy, N. Gisin, Phys. World (March) 

11, 41 (1998); E. Knill, R. Laflamme, G. Milburn, Na- 
ture (London) 409, 46 (2001). 
[3] P. Zoller et al, Eur. Phys. J. D 36, 203 (2005). 
[4] J. Wenger, R. Tualle-Brouri, P. Grangier, Phys. Rev. 

Lett. 92, 153601 (2004). 
[5] M. Hillery, R. F. O'Connell, M. O. Scully, E. P. Wigner, 

Phy. Rep. 106, 121 (1984). 
[6] H. J. Kimble, M. Dagenais, L. Mandel, Phys. Rev. Lett. 

39, 691 (1977). 
[7] R. Short, L. Mandel, Phys. Rev. Lett. 51, 384 (1983). 
[8] V. V. Dodonov, J. Opt. B: Quantum Semiclass. Opt. 4, 

R1-R33 (2002). 
[9] G. S. Agarwal, K. Tara, Phys. Rev. A 43, 492 (1991). 
[10] A. Zavatta, S. Viciani, M. Bellini, Science 306, 660 

(2004); A. Zavatta, S. Viciani, M. Bellini, Phys. Rev. 

A 72, 023820 (2005). 
[11] M. S. Kim, E. Park, P. L. Knight, H. Jeong, Phys. Rev. 

A 71, 043805 (2005). 
[12] C. T. Lee, Phys. Rev. A 52, 3374 (1995). 
[13] G. N. Jones, J. Haight, C. T. Lee, Quantum Semiclass. 

Opt. 9, 411 (1997). 
[14] A. I. Lvovsky, H. Hansen, T. Aichele, O. Benson, J. 

Mlynek, S. Schiller, Phys. Rev. Lett. 87, 050402 (2001); 

H. Hansen, T. Aichele, C. Hettich, P. Lodahl, A. I. 

Lvovsky, J. Mlynek, S. Schiller, Opt. Lett. 26, 1714 

(2001) ; A. I. Lvovsky, J. H. Shapiro, Phys. Rev. A 65, 
033830 (2002); A. Zavatta, M. Bellini, P. L. Ramazza, 
F. Marin, F. T. Arecchi, J. Opt. Soc. Am. B 19, 1189 

(2002) ; A. Zavatta, S. Viciani, M. Bellini, Phys. Rev. 
A 70, 053821 (2004); S. A. Babichev, B. Brezger, A. I. 



Lvovsky, Phys. Rev. Lett. 92, 047903 (2004). 
[15] J. K. Asboth, J. Calsamigila, H. Ritsch, Phys. Rev. Lett. 

94, 173602 (2005). 
[16] M. G. A. Paris, Phys. Rev. A 59, 1615 (1999); P. van 

Loock, S. L. Braunstein, Phys. Rev. Lett. 84, 3482 

(2000); S. Scheel, L. Knoll, T. Opatrny, D.-G.Welsch, 

Phys. Rev. A 62, 043803 (2000); M. S. Kim, W. Son, 

V. Buzek, P. L. Knight, ibid. 65, 032323 (2002); X. B. 

Wang, ibid. 66, 024303 (2002); 66, 064304 (2002); M. 

M. Wolf, J. Eisert, M. B. Plenio, Phys. Rev. Lett. 90, 

047904 (2003); J. S. Ivan, N. M ukunda, R. Simon, e- 

print arXiv : quant -ph/0603255 
[17] E. C. G. Sudarshan, Phys. Rev. Lett. 10, 277 (1963); R. 

J. Glauber, Phys. Rev. 131, 2766 (1963). 
[18] W. Schleich, J. A. Wheeler, Nature (London) 326, 574 

(1987); S. Schiller et al., Phys. Rev. Lett. 77, 2933 

(1996). 

[19] M. Hillery, Phys. Rev. A 35, 725 (1987). 

[20] C. T. Lee, Phys. Rev. A 44, R2775 (1991). 

[21] A. Wiinsche, V. V. Dodonov, O. V. Man'ko, V. I. 
Man'ko, Fortschr. Phys. 49, 1117 (2001); P. Marian, 
T. A. Marian, H. Scutaru, Phys. Rev. Lett. 88, 153601 
(2002); V. V. Dodonov, O. V. Man'ko, V. I. Man'ko, A. 
Wiinsche, J. Mod. Opt. 47, 633 (2000); P. Marian, T. 
A. Marian, H. Scutaru, Phys. Rev. A 69, 022104 (2004). 

[22] N. Liitkenhaus, S. M. Barnett, Phys. Rev. A 51, 3340 
(1995). 

[23] G. Vidal, R. F. Werner, Phys. Rev. A 65, 032314 (2002). 

[24] Trace norm of a partially transposed density operator 
g PT is given by ||£> PT ||i = 1 + 2N(p), where the nega- 
tivity N(p) is the sum [ Ai| of all the negative eigen- 
values of g PT . 



